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Abstract 

We investigate the adsorption of chlorine on the Cu(lll) surface with full 
potential all-electron density functional calculations. Chlorine adsorption at 
the fee hollow sites is slightly preferred over that at the hep hollow. The 
adsorption geometry is in excellent agreement with electron diffraction and 
ion scattering data. Adsorption energies and surface diffusion barriers are 
close to those deduced from experiment. 
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I. INTRODUCTION 



The determination of surface structure and adsorption energetics is a key issue in surface 
scienceS. Adsorption on metallic surfaces is of particular importance due to its relevance to 
many industrial processes. First principles simulation, based on density functional theory, 
is now one of the key tools for studying and developing an understanding of these systemsi. 

The adsorption of halides on metal surfaces and especially CI on Cu(lll) has been the 
subject of extensive experimental studies. The first study in 1977 reported on results from 
low energy electron diffraction (LEED), Auger electron spectroscopy (AES), work function 
measurements and desorption experiment si. Further work has been performed using surface 
extended X-ray absorption fine structure (SEXAFS) and photoelectron diffractioniH, normal 
incident X-ray standing wavefield absorption (NIXSW)&II, secondary ion mass spectrome- 
try (SIMS)0, scanning tunneling microscopy (STM)§, LEED, AES, thermal and electron 
stimulated desorptionEl. However, we are not aware of any first principles simulations of 
this system. This is no doubt in part due to the large computational effort needed to study 
the complex surface geometry with sufficient accuracy. In this article we present the results 
of an extensive study of the surface adsorption geometry and energetics within a gradient- 
corrected density functional formalism. It turns out that a very high accuracy is necessary 
to distinguish the different adsorption sites and thus careful calculations capable of resolving 
energy differences of the order of only a few meV are required. In view of this we present 
details of the tests performed to establish the numerical accuracy of the calculations. 

The outline of the article is as follows. In section |T|, we give the details of the computa- 



tional method. In section III, we present results for bulk Cu which document the accuracy of 



the numerical methodology and the functional used. In section |^ the surface geometry and 
energetics of the clean Cu surface are presented. The important numerical approximation 
of reciprocal lattice sampling is discussed in detail. In section we present a quantitative 
description of the adsorption of CI on the (111) surface of Cu. 
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II. METHOD 



The fundamental approximation made in density functional calculations is the choice of 
the effective one-particle potential. In this study we employed three distinct approxima- 
tions: the local density approximation (LDA) with Dirac-Slater exchangeEl and the Perdew- 
Zunger correlation functional^; the gradient corrected exchange and correlation functional 
of Perdew and Wang (GGA)li; and the Hartree-Fock (HF) approximation. 

A local basis scheme based on atom-centred Gaussian type orbitals was used. The 
calculations were performed with the CRYSTAL98 software^. The Cu basis set was chosen 
as [6s5p2(i]lll. The values of the inner exponents (a set consisting of one s contraction, 
3 sp contractions and 1 d contraction, i.e. [4s3]9lc/]), is given in table |. The outermost 
exponents were optimized in calculations on copper bulk at the GGA, LDA, and HF level. 
It is interesting to note that for Cu, in contrast to lithiumlli, it is possible to optimize the 
total energy with respect to the exponents without numerical difficulties even in the HF 
approximation. We believe that this is due to the more localized nature of the orbitals in 
Cu. 

For CI the [5s4p] basis set developed in previous calculations on alkali halides was used^, 
with outermost sp exponents of 0.294 and 0.090. An additional d-function with exponent 
0.5 was added yielding a [5s4pl(i] set. 

The density functional potential was fitted with an auxiliary basis set, which consisted, 
for both Cu and CI, of 12 s and p functions with exponents taken to be a geometrical 
sequence from 0.1 to 2000, and 5 d functions with exponents in the range 0.8 to 100. 

As in a recent study of lithiumil, we employ finite temperature density functional 
theoryS to ease the numerical integration over fc-space with the occupation calculated ac- 
cording to the Fermi function at finite temperature T (Ref. [2^ ) and the zero temperature 
energy is approximated byil ^{E{T) + F{T)); with E the energy and F the free energy 
F = E — TS, where S is the electronic free entropy. The important approximation associ- 
ated with the choice of the sampling nets employed in reciprocal space integration will be 
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discussed in detail below. 



III. BULK PROPERTIES 

We calculated the band structure and cohesive properties of Cu bulk in order to document 
the performance of our numerical approximations and to compare different choices for the 
effective one-particle potential (LDA, GGA, and HF). 

Figure displays the LDA band structure for bulk copper which is in excellent agreement 
with that reported in the literatureil. The LDA band structure remains virtually unchanged 
when only LDA exchange is included. Using GGA exchange and correlation leads to a similar 
band structure. The Hartree-Fock band structure (figure 0) reveals the expected pathology 
associated with the non-local exchange interaction and neglect of correlation in a metal. The 
overall bandwidth is too large by a factor of nearly two and the ci-bands too low in energy 
compared to the wide Asp band. 

The structure and cohesive properties are reported in table |I[ The GGA cohesive energy, 
lattice constant and bulk modulus are in good agreement with those observed. The LDA 
tends to over-bind the system resulting in a somewhat too low lattice constant and thus a 
high bulk modulus. The HF solution is dramatically underbound, the lattice constant far 
too large and thus the bulk modulus is very low. It is clear that the structure and energetics 
of the bulk crystal are best described within the GGA. 



IV. CU SURFACES 

We studied the clean Cu surface in some detail in order to establish the accuracy of 
predicted surface structures and energies. 

We used two methods to compute the surface energyili^l. Firstly, by extrapolating from 
calculations on slabs with a different number of layers {n and m): 

1 Ti 
Esurface = ^(^siafe(n) - {E,iab{n) - E,iab{n - m))—) (l) 



and secondly by using an independent bulk energy: 

E surface = -;^{E siabi'n) — E^uik X n) (2) 

where all the quantities E surface, Esiab{n), and Ebuik are expressed as energies per atom. 
In Ref. we gave a detailed comparison of the application of both formulas to the Li surface. 
The electronic structure of Cu is more complex than Li so as a first step we investigated the 
dependence of the surface energies on the electronic temperature and number of sampling 
points. 

We focus the discussion now on the Cu (111) surface. In figure we display the temper- 
ature and sampling point dependence of the surface energy obtained from equation |I] using 
two slabs with 3 and 4 layers. The density of reciprocal lattice points is determined by a 
shrinking factor. We used different shrinking factors of 2, 4, 8, 12, and 16 which result in 
2, 4, 10, 19, and 30 points in the irreducible part of the Brillouin zone. As expected the 
surface energy converges with respect to fc-point sampling at higher temperatures but as the 
temperature is raised the extrapolation to the athermal limit becomes less accurate. 

For the very high accuracy required in the current study we selected a shrinking factor of 
16 and a temperature of ksT = O.OlEh (i?/i=27.2114 eV). This converges the total energy of 
the bulk to better than 0.1 mEh, and that of the three layer copper slab to better than 0.3 



mEfi- In table |ITI|, we show the variation of the computed surface energy with slab thickness 
obtained both with equation |l| and The convergence of the data based on equation |^ 
demonstrates the independent convergence of the bulk and slab energies with respect to k- 
space sampling. This is an important prerequisite for accurate studies of chemical processes 
at surfaces. In many surface studies bulk cells and /c-space sampling are chosen to eliminate 
systematic errors in the definition of the surface energy. We note that this is not sufficient 
to guarantee the accuracy of surface properties. 

Table ^ contains the computed surface energies of copper. As found in earlier studies on 
the homogeneous electron gasEl and lithiumS, the LDA calculation results in a higher surface 
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gy than the GGA. Previous LDA results@'0 are rather scattered but in reasonable 



agreement with those computed here. 

In addition, we studied the relaxation of the Cu (111) surface. For this purpose, the top 
layer of various slabs of different thickness was allowed to relax with the layer spacing within 
the slab fixed at the bulk value of The results in table |V| demonstrate that convergence 

is achieved for slabs with 4 or more layers. The computed inwards relaxation of 1.0 % is in 
very good agreement with recent experiments using medium-energy ion scattering (MEIS)ii 
where a contraction of 1.0±0.4% was measured for the first interlayer spacing; it also agrees 
with that deduced from LEED experiments where inwards relaxations of 0.3±l%ii and 
0.7±0.5%i have been reported. A previous LDA calculation reported a relaxation of 

In this calculation the second interlayer spacing, which we kept fixed at the 
bulk value, was also reported to decrease by 0.64 % while an experimental studyil found a 
decrease of 0.2 %. We find that the relaxation of the top layer reduces the surface energy 
very slightly (by 0.003 ^). 

V. CHLORINE ADSORPTION ON THE CU(lll) SURFACE 

The adsorption of chlorine on the Cu(lll) surface was first studied by Goddard and 
Lamberti. At a coverage of one third of a monolayer, a structure displaying long range 
order in a VSx VS R30° pattern was observed and a desorption energy of 236 kJ/mol was 
obtained (see also Ref. 0). More recently the x ^/3 R30° pattern has been confirmed 
in STM imageJiil. The distance between the CI layer and the surface Cu layer has also 
been deduced from SEXAFS and photoelectron diffraction experimentsi'i. In this study, a 
distance of 2.39 ± 0.02 A between nearest Cl-Cu neighbours was obtained (corresponding 
to an interlayer spacing of 1.88 ± 0.03 A ); the adsorption site was identified as the fee (face 
centered cubic) hollow (see figure A distance of 1.81 ±0.05 A between the CI layer 

and a notional Cu-layer corresponding to unrelaxed bulk termination was deduced from 
NIXSW datal^i. Comparing the latter two results (SEXAFS and NIXSW) gives evidence 
for a slight inwards relaxation of the surface Cu layer similar to that observed on the clean 
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surface. SIMS measurementsc3 yield an interlayer spacing of 1.87 ± 0.04 A . Recently, in 
NIXSW measurements, the fee site was confirmed as the adsorption site, but also a small 
occupation of the hep (hexagonal closed packed) hollow was observed!. 

We optimized the structure of a slab with a chlorine coverage of one third of a monolayer 
in a X R30° cell. GGA calculations with the computational parameters described 
in section |^ were performed which led to 73 sampling points in the irreducible part of the 
reciprocal lattice. Slabs consisting of 3 and 4 layers of Cu were used in which the Cu atoms 
in the top layer and the CI layer were allowed to relax perpendicular to the surface. CI 
was adsorbed in a number of sites; the threefold hep hollow, the threefold fee hollow, a top 
position, and a bridge position (figure ||). 

The relaxed geometries and adsorption energies are presented in table 0. Adsorption 
on the fee or hep sites is clearly preferred over adsorption on bridge or top sites (by 3 and 
17 mEh, respectively). The Cl-Cu interlayer distance is 1.89 A for the fee hollow (1.90 
A for the hep hollow) in excellent agreement with that deduced from experiment. The 
next neighbour Cl-Cu distance is 2.40 A for the fee site, 2.41 A for the hep site, 2.33 A 
for the bridge site and 2.17 A for the top site. This is consistent with the idea that the 
next-neighbour bond strength is greater, and therefore the bond shorter, when fewer nearest 
neighbours are available!!. The inwards relaxation of the copper layer of ~ 0.04 A (1.9 %) 
is insensitive to the CI adsorption site and similar to that found for the clean surface (1.7 
%, table 0). 

The energy difference between the fee and hep sites is 0.3 mEh and increases to 0.4 
mEh when the Cu surface is relaxed. To investigate the dependence of this rather delicate 
result on the number of layers of the slab, we performed similar calculations for fee and hep 
site on a four layer slab. The fee site remains the preferred one by 0.2 mE^. The relaxed 
surface geometry (i.e. the interlayer distance Cl-Cu) is identical to that computed for the 
three layer slab. The relaxation of the top Cu layer reduces slightly to 0.03 A(1.4 %) which 
is a similar trend as for the clean surface (table 0). This value is within the errorbars of 
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the experimentally expected value. The optimization indicates that the fee site is becoming 
more stable relative to the hep site as the interlayer spacing between CI and the second Cu 
layer is reduced. This is demonstrated in Figure |^ where, at a fixed interlayer spacing of 
1.886 A between CI and the top Cu layer, the top Cu layer is allowed to relax inwards. 
As the relaxation increases, the fee site becomes more stable with respect to the hep site. 
The small energy difference between the fee and hep adsorption sites is in agreement with 
observations of predominantly fee-site adsorption accompanied by a small occupancy of the 
hep sitei. 

The calculated adsorption energy of 0.135 Eh (per CI atom) is in reasonable agreement 
with that deduced from desorption measurements (~ 250 kJ/mol corresponding to 0.095 
The experimental desorption energy is obtained from an Arrhenius model and 
therefore depends on a pre-exponential factor which was chosen to be lO^^i (^j^gf_ |||)^ but 
questioned in Ref. |^. 

The activation energy for surface diffusion was measured experimentally by electron 
stimulated diffusion0, i.e. the activation energy for CI atoms to diffuse into an area where 
CI depletion is induced by an electron beam. A value of 19 kJ/mol or 7 mEh was measured. 
A crude estimate from our calculations is possible if we assume that the energy difference 
of 3 mEh between the hollow site and the bridge site is of the order of the activation energy 
for diffusion. 



VI. CONCLUSION 

We have demonstrated that a quantitative description of the adsorption of CI on the 
Cu(lll) surface can be achieved with full potential, all-electron, GGA calculations. The 
surface structure and adsorption energies are in excellent agreement with experiment. The 
Cl-Cu bond length is found to be 1.89 A and the top Cu layer relaxes inwards by 1.4 %. The 
fee hollow is found to be the preferred adsorption site, with the hep site being ~ 0.2 mEh 
(5 meV) higher in energy. The energy difference between the two sites is sufficiently small 
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that it may alter slightly if a larger slab geometry was considered. Similarly, by analogy to 
the results for the clean Cu slab we might expect the inwards relaxation of 1.4 % to become 
slightly lower with larger slab geometries. The adsorption energy and surface diffusion energy 
is in reasonable agreement with those which can be estimated from experiment. Bulk and 
surface properties of Cu metal are found to be in very good agreement with experiment and 
previous calculations, when the gradient corrected GGA functional is used. 
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TABLES 

TABLE I. The [6s5p2d] copper basis set with the two outermost sp exponents and the outer- 
most d exponent optimized within the GGA. Those based on LDA and HF are shown in brackets. 



exponent 

398000.0 

56670.0 

12010.0 

3139.0 

947.2 

327.68 

128.39 

53.63 



s contraction 
0.000227 
0.001929 
0.01114 
0.05013 
0.17031 

0.3693 

0.4030 

0.1437 



p contraction 



d contraction 



sp 



1022.0 
238.9 
80.00 
31.86 
13.33 
4.442 



-0.00487 

-0.0674 

-0.1242 

0.2466 

0.672 

0.289 



0.00850 
0.06063 
0.2118 
0.3907 
0.3964 
0.261 



sp 



54.7 
23.26 
9.92 
4.013 



0.0119 
-0.146 
-0.750 
1.031 



-0.0288 
-0.0741 
0.182 
1.280 



sp 



1.582 



1.0 



1.0 
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sp 0.596 (LDA: 0.610; HF: 0.555) 1.0 1.0 

sp 0.150 (LDA: 0.150; HF: 0.170) 1.0 1.0 

d 48.54 0.031 

13.55 0.162 

4.52 0.378 

1.47 0.459 

d 0.392 (LDA: 0.392; HF: 0.423) 1.0 



14 



TABLE II. The ground state properties of bulk copper. Energies are in Hartree units, lattice 
constants in A, bulk moduli in GPa. 



a Ecoh B 

LDA 3.53 0.182 195 

GGA 3.63 0.143 155 

HE 3.95 0.018 69 

Lit. (KKR)il 3.59 0.155 158 

Lit. (LDA)il 3.62 0.133 147 

exp. 3.6040 0.12gil 1420 
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TABLE III. The convergence of the surface energy of the Cu(lll) surface as a function of the 
number of layers, computed using a shrinking factor of 16 and a smearing temperature of 0.01 Eh. 
For the Cu(lll) surface, 1 corresponds to 76.4 ^ for a lattice constant of 3.63 A. 

number Esiab{n) Esiab{n) - Esiab{n - I) E^^r/oce using E^^r /ace using 

of layers Esiab{n) - Egiabin - 1) Ebuik = -1640.698861 





Eh 


Eh 
atom 


atom 


Eh 
atom 


1 


-1640.653606 


- 


- 


0.02263 


2 


-3281.350305 


-1640.69670 


0.02155 


0.02371 


3 


-4922.049528 


-1640.69922 


0.02407 


0.02353 


4 


-6562.748468 


-1640.69894 


0.02365 


0.02349 


5 


-8203.447228 


-1640.69876 


0.02329 


0.02354 


6 


-9844.146061 


-1640.69883 


0.02347 


0.02355 


7 


-11484.844913 


-1640.69885 


0.02352 


0.02356 


8 


-13125.543718 


-1640.69881 


0.02336 


0.02359 


9 


-14766.242595 


-1640.69888 


0.02365 


0.02358 


10 


-16406.941437 


-1640.69884 


0.02349 


0.02359 


11 


-18047.640289 


-1640.69885 


0.02354 


0.02359 
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TABLE IV. The surface energy (^) of the low index copper surfaces. 

surface LDA GGA Ref. 1^ (LDA) ; Ref. || (LDA) 

(100) 2.52 2.01 2.09 ; 1.712 

(110) 2.67 2.15 2.31 ; 1.846 

(111) 2.26 1.80 1.96 ; 1.585 
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TABLE V. Relaxation of Cu (111) surface, as a function of the number of layers, computed 
within the GGA. The top layer is allowed to relax, the other layers are kept at a fixed distance of 
^'^^ corresponding to the bulk lattice constant 

number of layers relaxation in A (in %) 

3 -0.035 (-1.7 %) 

4 -0.025 (-1.2 %) 

5 -0.023 (-1.1 %) 



6 -0.022 (-1.0 



'0 



Literature: 

exp.i -1.0 ± 0.4 % 

exp.ii -0.3 ± 1 % 

exp.§ -0.7 ± 0.5 % 

LDAii -1.27 % 
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TABLE VI. Adsorption of CI on the Cu(lll) surface. S1-2 is the change in interlayer spac- 
ing between first and second copper layers (in A) relative to the bulk value, dci-Cu top layer is 
the distance between the CI and top Cu layer (in A). The adsorption energy is the difference 

Ec\ at Cu(lll) - -Scu(lll) - Eci- 



fee site 
n 
3 
3 
4 

hep site 
3 
3 
4 

bridge site 
3 

top site 
3 



Sl-2 

(no relaxation allowed) 
-0.043 
-0.032 

(no relaxation allowed) 
-0.038 
-0.031 

-0.040 

-0.041 



C?Cl-Cu top layer 
1.89 
1.89 

1.89 

1.90 
1.90 
1.90 

1.94 

2.17 



Eadsorption in Eh, per CI atom 
0.13542 
0.13566 
0.13583 

0.13510 
0.13527 
0.13563 

0.13265 

0.11888 
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FIGURES 



FIG. 1. LDA band structure at the equilibrium lattice constant. 

FIG. 2. HF band structure at the equilibrium lattice constant. 

FIG. 3. (Ill) Copper surface energy with four different reciprocal lattice samplings extracted 
from two slabs with 3 and 4 layers using ^{E{T) + F(T)) and equation ^ 

FIG. 4. The structures considered for CI, adsorbed on the Cu(lll) surface at one third coverage, 
in a \/3 x \/3 R30° unit cell. When CI is adsorbed in an fee hollow, it sits above a Cu atom in 
the third layer (upper left) while in the hep hollow it is above an atom in the second layer (upper 
right). The top position is CI adsorbed vertically above a surface atom (lower left). In the bridge 
position it is vertically above the middle of two surface atoms (lower right). 

FIG. 5. Total energy and energy difference of fee and hep adsorption site as a function of the 
inwards relaxation of the top Cu layer. The data in this figure is from a four-layer slab at a fixed 
interlayer spacing between CI and top Cu layer of 1.886 A for both fee and hep adsorption site. 
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Fig. 1 
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Fig. 2 
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Fig. 4 
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CI on Cu(111) surface 
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Fig. 5 
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